A microscopic theory for the incommensurate transition in TiOCl 
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We propose a microscopic mechanism for the incommensurate phase in TiOX compounds. The 
model includes the antiferromagnetic chains of Ti ions immersed in the phonon bath of the bi- 
layer structure. Making use of the Cross-Fisher theory, we show that the geometrically frustrated 
character of the lattice is responsible for the structural instability which leads the chains to an 
incommensurate phase without an applied magnetic field. In the case of TiOCl, we show that our 
model is consistent with the measured phonon frequencies at T = SOOiC and the value of the incom- 
mensuration vector at the transition temperature. Moreover, we find that the dynamical structure 
factor shows a progressive softening of an incommensurate phonon near the zone boundary as the 
temperature decreases. This softening is accompanied by a broadening of the peak which gets asym- 
metrical as well when going towards the transition temperature. These features are in agreement 
with the experimental inelastic X-ray measurements. 

PACS numbers: 63.20.kk, 75.10.Pq, 75.50.Ee 
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I. INTRODUCTION 

In 1955, when writing an introductory solid-state 
textbook,^ Rudolph Peierls discovered that no one- 
dimensional metal could exist due to the electron-phonon 
coupling. Indeed, a half-filled metal is unstable towards 
a lattice dimerization. The system undergoes a metal- 
insulator transition because the loss of elastic energy is 
made up for by the electronic energy acquired when the 
gap is opened. A quite similar process takes place in a 
one-dimensional antiferromagnet, giving rise to the so- 
called 'spin-Peierls' (SP) transition.^ The lattice dimer- 
izes and the antiferromagnetic quasi long range order 
is replaced by a gapped singlet state. From this point 
of view, no one-dimensional antiferromagnetic system 
should be stable. 

In spite of this general prediction, only few spin-Peierls 
materials were found. Most of them were observed in the 
early 70s in some organic charge transfer systems. How- 
ever, direct spectral characterization was lacking due to 
the unavailability of large crystals for inelastic neutron 
scattering measurements. In this context the discovery 
of CuGeOa in 1993, the first inorganic spin-Peierls sys- 
tem, renewed the interest in this subject. Large enough 
crystalline samples of CuGeOa could be obtained to un- 
dertake a detailed spectral characterization. An impor- 
tant conclusion of these studies was that in CuGeOs, the 
nonadiabatic character of the phonons^ is so important 
that the SP transition is not driven by a softening of 
a precursive phonon mode. It has been shown that an 
extension of the canonical theory of Cross-Fisher^ could 
explain this feature.^ Moreover, the dispersion relation of 
the phonons in the direction perpendicular to the mag- 
netic chains should be taken into account in a nonadia- 
batic SP system like CuCeOs-^ 

The recent discovery of a spin-Peierls transition in the 
TiOX (X = CI, Br) family has opened new questions 



about this spin-Peierls paradigm. The essential build- 
ing blocks of these compounds are bilayer structures of 
magnetically active Ti ions connected by O ones. The 
position of a Ti ion in a layer is shifted with respect to 
the other in the neighboring layer, forming something 
like an anisotropic triangular structure. In fact TiOCl 
was initially thought to be a candidate for a resonating 
valence bond state, ^ but it has been found to be mainly 
a one-dimensional magnetic system^ with the Ti dxy or- 
bitals pointing toward each other in the crystallographic h 
direction.^ The high temperature magnetic susceptibility 
is well described by the Bonner-Fisher curve, indicating 
a nearest neighbor magnetic exchange J ~ 660i^.^ 

The phase diagram of TiOCl does not correspond to a 
canonical SP system in the sense that an incommensurate 
intermediate phase appears between the high tempera- 
ture uniform phase and the low T dimerized phase j^ii^ 
The transition temperatures are found to be Td ~ 66i^ 
and Tc2 ~ 92i^. In Ref. 10 it is also reported a very 
large energy gap of about A?>{)K in the low temperature 
phase and a pseudo spin gap below 135i^r. The order of 
both transitions is still under debate J^ii^ii^ 

The origin of the intermediate phase is controversial 
and not yet well understood. As the position of a Ti 
ion in a chain is shifted with respect to another one in 
the neighboring chain, it is plausible to speculate that 
some type of competition between the in-chain and out- 
of-chain interactions could be the origin of the incom- 
mensurate phase. This was the idea of Riickamp et al}^< 
who proposed a Landau theory for the incommensurate 
transition. This phenomenological theory includes the 
tendency to dimerize of each chain and the coupling of 
the order parameter with the neighboring chain. As a 
result the incommensurate phase is accounted. However 
the sign and the value of the parameters are phenomeno- 
logically chosen and the connection with the underlying 
microscopic theory is not clear. 
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In view of the previous discussion, in this paper we 
propose a simphfied microscopic model which accounts 
for the transition from the uniform to the incommensu- 
rate phase in TiOCl. Our model contains the relevant 
magnetic interaction, the phonons and the spin-phonon 
coupling. The paper is organized as follows: In Sec. [Ill 
we present the model and show that the essential ingre- 
dient leading to the incommensurate transition is that 
the phononic dispersion of the modes near tt is linear. In 
Sec. Uniwe apply the model to TiOCl. We start by fit- 
ting the parameters in order to account for the measured 
phononic frequencies and the value of the incommensu- 
rate wave vector at the transition temperature. After 
that, we obtain the dynamical structure factor of the 
phonons and compare it with the one obtained in X-ray 
experiments. The softening of a phonon and the loss of 
spectral weight when the transition approaches is cor- 
rectly accounted. Finally, in Sec. IIVI we discuss our 
results in comparison with another approach and sum- 
marize them. 



II. MODEL AND INCOMMENSURATE 
TRANSITION 

We include the Ti atoms on the bilayer structure which 
interact by harmonic forces as depicted in Fig. [H For 
simplicity we consider ionic displacements only in the 
direction of the magnetic chains. As the measured mag- 
netic susceptibility of TiOCl is well reproduced by a ID 
Heisenberg model, we assume that the Ti atoms are mag- 
netically coupled only along the b direction. Moreover, 
as a direct exchange seems to be the dominant Ti-Ti in- 
teraction, we assume that J{Au) is modulated by the 
movement of nearest neighbor Ti ions in the chain direc- 
tion. Our spin-phonon Hamiltonian reads: 



H — Hph -\- Hs -\- Hgph 



(1) 



P2. 



^ inter 



[( 



Hsph = ^<^('^i+l,j — Uij) Sij • Si+ij, 



where Pij is the momentum of the atom i of the chain 
j, Uij are the displacements from the equilibrium posi- 
tions along the direction b of the magnetic chains, Sij are 
spin- 1 operators with exchange constant J = J{Au = 0) 
along the b-axis of a non-deformed underlying lattice, 
a = {dJ{Au)/dAu)\/\u=o^ and Kin and Kilter are the in- 
chain and inter-chain harmonic force constants as shown 
in the Fig. [H 



J(Au) 




FIG. 1: Schematic representation of our simplified model. 
Only Ti atoms are included over the ab plane. Kin and Kinter 
are the harmonic force constants acting when two neighbor 
atoms in the same chain and in neighboring chains respec- 
tively, move from their equilibrium position. 



Representing the atom displacements and momenta 
through phonon normal coordinates, we can write the 
phonon dependent parts of the Hamiltonian as: 

q ^ ^ 

where ^^(q) is the phonon dispersion relation which reads 



^'(q) = -^[K^n sin' ^ + K,nter{l - COS ^ COS ^)], (2) 

for our model and g = The quantized phonon nor- 

mal coordinates are defined by 



Q(q) = 



(a-q + aq). 



It will be also useful to introduce a dimensionless spin- 
phonon coupling constant defined as in Ref. i: 



TTjQ?{qy — tt) * 



To compare with X-ray scattering data, the dynami- 
cal structure factor 6'(q, uj) has to be determined. It is 
related to the phononic retarded Green's function^^ by 



^(q,c^) 



l](q) Im P^^^(q,cj) 
1 - e^^ ' 
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where 



DO(q,cj + i^) 



(3) 



1 - DO(q, uj + i(^)n(q, uj^id)' 

denotes the non-interacting phonon Green's function 
defined by 

n 2 
V (q,icJn) = J. 7^ F^^TT' 

and n(q, cc;) = 11(^2/, cj) is the phonon self-energy. In 
the following, we treat the spin-phonon interaction using 
a random phase approximation (RPA) for the phonon 
self-energy term along with the expression obtained by 
Cross and Fisher^ using bosonization for the dimer-dimer 
correlation function in the Heisenberg model: 



with 



0.37 



T 



1(1 



2^ 



X h 



uj - ^Jb{qy - tt) 



h{k) = 



^4 



87rr(| + ^i/c) 



Let us start by identifying the structural transition of 
the model by searching the poles of the retarded Green's 
function when uj = over the (qx^Qy) plane of the first 
Brillouin zone. They are given by the equation of the 
zeros in the denominator of Eq. ([3|) which reads: 



n\q)^2U{qy,0) = 0. 



(4) 



From this equation we obtain T(q), the temperature 
where the renormalized frequency vanishes for each q. 
The highest of these temperatures signals the q wave 
vector of the structure that will be developed at that 
T. For an isolated chain this transition takes place at 
q = (0,7r), which is the usual spin-Peierls transition to- 
wards the dimerized phase. 

How does this situation change when the system ac- 
quires a transversal dispersion due to an elastic coupling 
between the chains? Certainly, if the chains arrange in a 
rectangular geometry (atoms in different chains align in 
phase), the dispersion of the phonons in the chain direc- 
tion has a local maximum at = tt and the T(q) given 
by dH) has a maximum at q = (0, tt). Thus, no complete 
softening of an incommensurate mode is found. 

A crucial consideration here is that the dispersion given 
by Eq. ([2]), when expanded near % = tt, reads: 



M 
TT-qy 



[Kin H~ Kinter 



o cos— )], 



i.e. it increases linearly near tt. Plugging this expression 
in Eq. (|4]) and expanding up to order 5 we obtain: 



T{5) 



g'^M\ao\ 



^ ^ 1 Kinter COs(^)(5 
2 Kinter ~t~ Kin 



(5) 



where it was enough to calculate n(g^,0) at qy = ir as 

n(7r,0) = with ao = -0.257743, because the first 

correction is of second order in S. 

The temperature given in ([5]) decreases with q^ and 
increases for increasing 5. Therefore an instability is pre- 
dicted dit qx = and qy shifted from tt. The next orders 
of T{S) stabilize a local maximum at a given S. Thus, 
due to the linear dispersion of the relevant phonon mode, 
the model ([T]) has a transition from a uniform to an in- 
commensurate phase by lowering the temperature. The 
physical reason why this transition appears in the trian- 
gular arrangement and not in the rectangular one is that 
the elastic energy that the system gains by distorting 
the lattice in a modulated pattern decreases faster in the 
first case, when the wave vector of the modulation moves 
away from tt. This reduction overcomes the cost in the 
magnetic free energy the system should pay to separate 
from the dimerized state. This reduction does not take 
place in the rectangular case and the system never goes 
to an incommensurate phase but directly to a dimerized 
one. 



III. APPLICATION TO TIOCL 

Having shown that the uniform-incommensurate tran- 
sition is present in our model, we will try to fit the pa- 
rameters to account for the available experimental data of 
TiOCl. We intend to reproduce the experimental incom- 
mensurate wave vector qy 3.04^^^^^ and the renormal- 
ized frequency of the phonons near the zone boundary at 
300i^ obtained in Ref. 15 by inelastic X-ray scattering. 
In order to do that we have to fix the free parameters of 
our model, i.e. the elastic constants and the spin-phonon 
coupling. Usually the critical temperature is used to fix 
the spin-phonon coupling. This value could be obtained 
by numerical solution of Eq. [4] once the other parame- 
ters are known. The values of Kin and Kinter are fixed 
in such a way to fit the phononic frequencies measured 
by inelastic X-ray scattering at T = 300i^ with the ones 
obtained from the positions of the peaks of S{q^u;). 

The va lue of the in-chain elastic constant was set to 
Qin = y^4:Kin/M = 6.6 meV and the inter-chain one 
was fixed to Kinter / Kin = 3.38. This enables us to fit the 
experimental frequencies near qy = tt dit 300K with an 
spin-phonon coupling constant A = 0.17. The transition 
temperature is at 92K as in TiOCl. However, using these 
parameters we obtain an incommensurated wave vector 
qinc = 3.125, i.e. it is close to tt and far away from 
the experimental value. Our results together with the 
experimental ones are shown in Fig. O In this figure 
we see that, as we move away from the zone boundary, 
the dispersion curve quickly decreases and separates from 
the experimental points. Moreover, the bare dispersion 
curve is very close to the 300i^ one. 

We have obtained a quite small value of the incom- 
mensuration with those parameters. Therefore, in or- 
der to reproduce the experimental one, we increase the 
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FIG. 2: (Color online) Evolution with the temperature of 
the dressed frequency of the phonons along the (0,^^) path. 
The model parameters are 4:Kin / M — 6.6 mel/, Kinter ~ 
3.38Kin and A = 0.17. The black dots are the experimental 
points obtained by inelastic X-ray scattering at T = 300i^.^^ 



spin-phonon coupling A. The new value turns out to be 
A = 0.41. This change has the undesired effect of in- 
creasing the transition temperature up to T^^c ~ 225i^, 
more than twice the experimental transition tempera- 
ture. Moreover, if we want to keep adjusting the ex- 
perimental X-ray phonons at 300 we have to increase 
the value of the bare phonon frequency because of the 
larger spin-phonon coupling. The chosen values are 

10.5 meV, Kinter 



kept the ratio K^r.+^^IK^. 
she □ 



3.38Kin, i-e. we have 
unchanged. The results are 




FIG. 3: (Color online) Evolution with the temperature of 
the dressed frequency of the phonons along the (0,^^) path. 
The curves were obtained from the positions of the peaks of 
th e dynam ical structure factor. The fitted parameters are 
yjAKin/M 10.5 meV, Kilter ^ 3.38Kin and A 0.41. 

curve does not fall rapidly as before and that the ex- 
perimental point which is the most separated from the 
zone boundary is also well adjusted. Nevertheless, we 
have noted that with this set of parameters the transi- 
tion temperature is overestimated. This problem is also 
present in the calculations of Abel et al}^ through an un- 
derestimation of the magnetic coupling J. We think that 



this disagreement is due to the adiabatic treatment of the 
phonons and will be solved by an extension of the results 
of Ref. 6 for the geometry of TiOCl. In this paper it was 
shown that a theory going beyond the mean field RPA 
treatment of the phonons by including nonadiabatic ef- 
fects predicts a reduction of the critical temperature for 
a given A. We plan to undertake this generalization in a 
forthcoming paper, but let us make a rough estimation of 
the renormalization of the transition temperature using 
formula (35) of Ref. 0. Note that the parameter flinter in 
Ref. is not exactly the same as in the present paper. It 
is in this sense that the following results should be taken 
as an estimation. Proceeding in this way we obtain: 



T = Tad{l 



inter / ^in)'^ 



117K, 



where Tad = 225K is the transition temperature obtained 
by the adiabatic approach. The non-adiabatic tempera- 
ture is now very close to the experimental one. Further- 
more, let us estimate the energy gap by combining Eqs. 
(26) and (32) of the same reference: 



3.04AJ 1 



-\- (Qjinter / ^in)^ 



i.e. we obtain an excellent agreement with the experi- 
mental gap of 430K.^^ From these estimations it can be 
inferred that the non-adiabatic treatment should solve 
the problem. 

For the present purpose, let us use the second set of 
parameters. In Fig. [H we show a plot of the dynam- 
ical structure factor along the (O,^'^) path of the first 
Brillouin zone for a range of frequencies and different 
temperatures. 

When the temperature is high [Fig. HJa)], we obtain 
peaks which reproduce the bare dispersion relation of the 
model. With decreasing temperatures [Fig. [H^b) to (d)], 
we observe the trace of the phonon softening along with 
a broadening and a height reduction of the peak. This 
is due to a redistribution of the spectral weight which 
is concentrated again in a defined peak at the transi- 
tion temperature for = and qy = 3.047 [Fig. \^d) 
and zoom in (e)]. This broadening of the peak indicates 
that phonons are overdamped by the interaction with the 
magnetic excitations. This behavior has been indeed ob- 
served experimentally in Ref. 15 [see Fig. 9(b) in this 
reference] and considered there as a disagreement with 
the Cross-Fisher theory. Instead, our results show that 
due to the interaction of the bare phonons with the two- 
spinon continuum of the magnetic subsystem, the dy- 
namical structure factor of the phonons does not show a 
Lorentzian character. It is indeed asymmetrical towards 
low frequencies (Fig. [5|). This behavior should be checked 
in future experiments. 
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FIG. 4: (Color online) Contour curves of the dynamical struc- 
ture factor as a function of qy {qx = 0) and u for different 
temperatures, (a) lOOOii^, (b) 300K and (c) (d) 225K 

(transition temperature) . In (e) we show a zoom of (d) around 
the transition point. Brighter colors correspond to higher in- 
tensities. We see that there is a softening of a group of modes 
near the zone boundary. In this region, there is a broaden- 
ing of the peak width with decreasing temperature due to a 
redistribution of the spectral weight. 



IV. DISCUSSION 

In this paper we have shown that a simple model 
including the elastic frustration of the bilayer structure 
is able to produce the incommensurate phase as seen 
in TiOCl. On the other hand, from an ab-initio DFT 
theory, it was recently proposed that a weak ferro- 
magnetic inter-chain coupling could be the origin of the 
incommensuration in this compound. As no quantitative 



estimation of the incommensurate wave vector has been 
done, we cannot compare the relative importance of 
each mechanism. More importantly, the actual values 
and the signs of the inter-chain exchanges remain as 
a controversial question and the results of the present 
work show that the elastic frustration in the interchain 
coupling is enough to produce the incommensurate 



phase. 




5 10 15 20 25 
fj (meV) 

FIG. 5: (Color online) Dynamical structure factor at qinc 
for different temperatures. At high T, we see a very definite 
peak which broadens and gets asymmetrical as temperature 
decreases. At T = 225K (the transition temperature), we 
observe the divergence of the peak for = 0. 



In summary, we have developed a model which explains 
the uniform-incommensurate transition in TiOCl. It is 
based on the competition between the tendency to dimer- 
ize of isolated antiferromagnetic chains and the frustrated 
elastic inter-chain coupling. The model predicts the soft- 
ening of an incommensurate mode. Moreover, it is in 
agreement with the loss of a coherent spectral weight for 
the phonons near the zone boundary as observed by in- 
elastic X-ray scattering. 
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